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' We consider the spectral function (/i > 0) for the Sturm-Liouville equation y + (A — q)y = 

(— I , on [0, oo) with the boundary condition y{0) = and where q has slow decay 0{x~") (a > 0) as x oo. 
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Abstract 
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2P ■ 1 Introduction 

. In a recent paper |]5| , we gave a new formula 



p (//) = TT ^scxp(— s ^ / q{x) sin 26 {x, ii)dx) (1. 1) 







X 

' for the derivative of the spectral function associated with the Sturm-Liouville equation 

y"{x) + {X-q{x)}y{x)^0 (0 < x < oo) (1.2) 
and the Dirichlet boundary condition 

y(0) = 0, (1.3) 
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the formula applying to the situation where 

q{x) e L(0,oo). (1. 4) 

In (|l. ij. > 0,s — ^JJi and Q[x^ [i) is the solution of the first-order differential equation 

Q {x,p.) — s — s^^qix^svY? Q{x,ii) (1. 5) 

such that 

0(O,m)=O. (1.6) 
Also, if ( |l. 3| ) is replaced by the usual general condition 

y(0) cosa + y (0) sina = (0 < a < tt), 

the initial factor s in ( |1. 1| ) is replaced by 

s csc^ a(s^ + cot^ 

and ( |1. 6| ) is replaced by 

6»(0,^) = -tan"^(stanQ;) (-tt < 6'(0,^) < 0). 

In this paper however, we keep to (|l. 3| ) purely for simplicity. 

Following § and (l8| (see also |]), we say that the problem ( |l. 3[ ) exhibits spectral concentra- 



tion at a point /io (> 0) if p has a local maximum at fiQ. Then p itself has a relatively sharp increase at 
fiQ. In lH we used ( |l. l\ ) to develop computational procedures, based on ( |l. 5| ) and ( |l. 6| ), for locating 
spectral concentration points jiQ with the emphasis on potentials q which decay rapidly as x —>■ oo. In 
particular, we identified a transitional property of 9(x,fi) as /i increases through fiQ which provides a 



sensitive test of even slight spectral concentration sections 2-3]. This property also features in this 
paper, and we give details of it later at the end of section 4. 

A typical example considered in ||, section 3.1] and suggested by j20[ example 166] is 

g(a;) = -ce-^/^cosx (c> 0) (1.7) 

with exponential decay. However, in [|[ section 5], we pointed out that our procedures are less reliable 
in cases of slower decay such as 

q{x) = -c(l + x)""cosa; (1. 8) 

with a > 1, and we also raised the further question of what can be said when < a < 1, in which case 
( p_. 4| ) fails and it is no longer clear that ( |l. l| ) is available. 

In this paper, we deal with these two outstanding matters. First, in section 2, we show that 
I 1\) continues to hold in certain situations where ( |1. 4D fails and q is only conditionally integrable on 
(0, oo). Certain values of have to be avoided because of the possibility of embedded eigenvalues and 
discontinuities in p but, in /i - intervals which avoid these values, the spectrum is absolutely continuous. 
Then, in section 3, we show how to extend the computational procedures in to locate reliably spectral 



concentration points for slow-decay examples such as ( 1. 8 ), with a < 1 allowed. 

Finally in this introduction, we note two other approaches to spectral concentration. The software 
package SLEDGE ]P,jl^, |2^ replaces q in ( |l . 2| ) by an approximating step-function q over a large 
interval (0,6). It then computes the spectral fimction p(q,b,^) for ( |l. 2| ) with q and the boundary 
conditions ( |l. 3| ) and y{b) ~ 0. When refinements of q and b lead to a stable output, the step- function 
p(g, 6, fi) provides an approximation to p{fi) to within a prescribed tolerance. SLEDGE is not restricted 
to ( |1. 4 ) but, when more explicit formulae such as ( |1. l\) and (1. 5) are available, these formulae provide 
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a more sensitive means of detecting spectral concentration 

The other approach to spectral concentration lies in the wider theoretical context of quantum reso- 
nances and spectral stability p^ . The resonances are non-real singular spectral points and, in a future 
paper, we intend to develop the connection between these non-real points and the real spectral concen- 
tration points fomid in this paper and in 



2 Absolute continuity of the spectrum 

The theory from which ( |l . l| ) is derived in Q was developed in the original work of Titchmarsh p3| , 
section 5.7] and Weyl p. 264]. This theory uses the fact that, subject to ( |1. 4| ), there are solutions 
of ( |1. 1[ ) which together with their derivatives are asymptotic to exp(±ia;\/A) and ±i\/Xexp(±ix\/X) as 
a; — > oo. More recently, the existence of solutions with these or similar asymptotic forms, irrespective of 

4[ ), provides an application of the subordination theory of Gilbert, Pearson and Stolz and 
leads to the absolute continuity of p{p,) in appropriate intervals. Such asymptotic forms are obtained by 
transforming ( |l. 2| ) into a first-order differential system to which the Levinson asymptotic theorem 
section 1.3] is applicable, and the necessary transformation methods were developed by Harris, Lutz and 
Eastham 10, In particular, Behncke 0], 0, |§] used the transformations in to establish 

the absolute continuity for potentials such as ( |l. 8| ) when a > 1/3 provided that certain resonance values 
of /i are avoided. In this section, we use the transformation due to Eastham and McLeod sections 
4.6-4.7], |l^ to both extend this result and establish ( |l. l\ ) for aU a > 0. 

In P, (4.1.8)] (|l. 2[) is considered with A = 1 and therefore some minor changes are required in the 
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transformation theory as presented in section 4.6]. As in j|, section 4.1], we take q to have the form 



q{x) = £.{x)p{x) 



(2. 1) 



where p has period 27r and 



£,{x) {x-^oo), ^ (x) e L{0, oo). 



Also, 



e(a;)^i*'(0,oo), e(^)eL*^+^(0,oo) 



(2. 2) 



for some integer M (> 1). Thus £,{x) — {1 + x) " (0 < a < 1) is the simplest example. The appropriate 
formulation of ( |l. 2| ) as a system is 



W (x) = {iAo + Rix)}W{x) 



(2. 3) 



where, as in ^, section 4.1], 



Ao = dg(VA, -VA), Rix) = ^i^{x)p{x)Dn 



(2. 4) 



with D = dg(i, — i) and fl having all entries unity. The connection between y and W is 



iVA — i\/A 



The resonance set a in (4.1.18)] is now replaced by 



(2. 5) 



a = {7VV4;iV=l,2,...}. 



(2. 6) 



We can now proceed as in j^. Lemma 4.6.1] with a transformation 



W = {exp(ePi + + ... + rPM)}Z, 



(2. 7) 



where the matrices Pm have period 2tt, dgPm = 0, and M is as in (2. 2). This takes (2. 3) into 



Z = {lAo + CAi + ... + r Am + S)Z, 



(2. 8) 



where S G L{0,oo) and A™ is diagonal with period 27r. Here A is excluded from the set a in (2. 6) 



Further, since R has trace zero in (2. 4), it follows that 



trA„, = -trF,„ = (m > 1) 



as in 1^, Lemma 4.9.2 (i)]. Thus A„i has the form 



A„ = dg(A('"),-A(™)). 



Finally, when A is real and positive, the conditions of [^, Theorem 4.6.1] are satisfied by ( 2. 3 ) and 
I 4|) , and then the A^ are all pure imaginary. 

The method used in section 5.7] for obtaining the Titchmarsh-Weyl m(A) function and the 

spectral function can be adapted to the situation which we have now in (|~1|)-(|7J). We note that 
( ^. 7| ) has the form 

W={I + Q)Z, (2.9) 
where Q{x) — o(l) {x oo), and we choose X so that {I + Q)^^ exists in [AT, oo). We also note that 



IS 



Z = (A + 



(2. 10) 



where A has the form 



A = dg(i^, -u) 



(2. 11) 



with 



i^^i\^ + o{l) (.T ^ oo). 



(2. 12) 



Further, v is pure imaginary when A is real and positive. 



The usual integral form of (2. 10) is 



Z{x) = $(a;)Z(X) + / ^{x)'^-^{:t)S{t)Z{t)dt 
Jx 



ith 



$(a;) = exp |dg i'(t)dt,- J v 



{t)dt 



(2. 13) 



Then, by ( |2. 9[ ), the solutions of ( |2. 3D satisfy 



where 



Wix) = {I + Qix)}<S>{x){I + Q{X)}-^W{X) 



+ {I + Q{x)} / ^{x)^-\t)S{t)W{t)dt, 



S = S{I + Q)-^ e L{X,oo). 



(2. 14) 



(2. 15) 



If A and \/A satisfy < argA < tt and < arg\/A < ^tt, it follows immediately from ( 2. 12 )-( 2. 15 ) 
and a Gronwall inequality that W^(a;)exp (J^ v{t)dt) is bounded on [X, oo). To use this property in the 



integral term in ( 2. 14 ), we define Ii — dg(l, 0) and I2 ~ dg(0, 1) in order to split the two entries in $ in 
( ^. 131 ). Then, by ( |2. 14D , we have 

W{x) = {hB{\) + o(l)}exp {- J v{t)dt^ {x 00) 

when A is non-real, where 

B{X)^{I + Q{X)}-'W{X) 

(of. 11^, (5.7.5)-(5.7.8)]). Also, when A has a real and positive value /i, ( ^. 14 ) again gives 

W{x) = /iA(^)exp ( / v{t)dt 



exp 
X \Jx 



iy{u)duj S{t)W{t)dt 



\Jx 



+ l2B{n)expi^- J iy{t)dtj+o{l), 

where v is now pure imaginary and A{fi) is the same as -B(/i) but with —ly instead of ly (cf. [p3| (5.7.2)- 
(5.7.3)]). Finally, in terms of the first component a(/i) of A(/i) and the second component 6(A) of B{X), 
the transformation ( 2. 5) ) back to y gives 



y(a;) = {6(A) + o(l)}exp ( -y iy{t)dt 



for A non-real and, when X = fj,, 



y{x) = a{fj.)exp / iy{t)dt + b{fi)exp - / iy{t)dt + o(l) 



X 



y (x) = i^/fla{ii)exp ( ^ iy{t)dt j — iy/JIb{fi)exp 



X 



v{t)dt +o(l). 



(2. 16) 



(2. 17) 



We now have the same type of asymptotic formulae as in |2^, section 5.7], from which (|l. l| ) follows as 
in Q. We indicate the details briefly, the only proviso being that /i ^ cr in ( 2. 6| ) as already mentioned. 
Let yi{x, A) and ^2(2;, A) be the solutions of ( |l. 2| ) which satisfy the initial conditions 

2/i(0,A) = l, 2/;(0,A) = 0, 2/2(0,A) = 0, y;(0,A) = l, 



and let ai, 61, 02, 62 denote the corresponding multipliers as in ( ^. 16 ) and ( ^. 17| ). Then yi my2 G i-^(0, 00) 
gives 

to(A) = -6i(A)/62(A) (imA ^ 0) (2. 18) 



as in |23, (5.7.9)]. When X = fj,, yi and 2/2 are real-valued and hence 
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in ( |2. 17D . Then W{yi,y2) = 1 gives 

im(Mi)(/i) - -(4Vm)"'- 



Hence, by (|2. 18|) , 
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lim im to(A) = {Ay/Jl \ 62(m) 1} 

X—^fi 

= lim ^^^-1/2^^2 

by ( |2. 17| ) and follows as in ||. 

Finally in this section, we note that the requirement fJ- ^ <y can be relaxed as follows when M = 1 in 
(|~2). Let c„ {—oo < n < oo) denote the complex Fourier coefficients of p{x) in ( |2. l| ). Then it is shown 



in 1^, section 4.2] ( see also [Q) that ( ^. 17] ) continuous to hold when /i = jN^ for some N provided that 
Cat = 0. Thus, altogether, ( |1. 1| ) is valid 

1. for all /i > when ^ e i(0,oo), 

2. for all /Lt > except those ^ = jN^ for which cjv 7^ 0, when M = 1 in ( |2. 2| ), 

3. for > and ^ cr when M > 2 in(|~J). 



3 An integration algorithm 

We aim to compute p in ( |l. l[ ) to within a reasonable degree of accuracy such as 10~^. The error in p is 
of course made up of a number of components: the truncation error due to the approximation of the semi 
infinite interval by a finite interval, the error inherent in the solving algorithm of the differential equation 
and the rounding error due to the rational approximation of real numbers in computer arithmetic. We 
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focus here on the truncation error and rely on the standard theory for both the error in the Runge Kutta 
algorithm as well as the floating-point numerical software on our computer system. The infinite integral 
m ( p~l| ) is of course truncated at a suitable value Xq. However, to achieve a truncation error of 10 in 
the integrand when, for example, a = 3 in (|l. 8| ) requires Xq — 10^, and integration over the large range 
(0, 10"^) is unreliable. The situation is much worse for smaller values of a. In this section, we develop 
an iterative algorithm which accelerates the convergence of the integral in ( |l. l|) and enables us to cope 
with potentials such as ( |l. 8| ) when a > 0. 

Guided by the example ( |l. 8 ), we give the algorithm for the case p(x) = cosx, so that (2. 1) is 



q{x) — £_{x) cos X (3. 1) 

and ^{x) is as before but with M > in (|~|). Our methods also cover the more general situation 
where p{x) is a finite Fourier series, but the details become more complicated. We require the following 
trigonometric identity, valid for any a, /?, 9 and x: 

8 cos X sin^ 9 sin(a6' + (3x) 
= 2 sin{a6i + (/3 + + 2 sm{a9 + (/3 - l)x} 

- sm{{a + 2)9 +{(3+ l)x} - sin{(a - 2)61 + (/? - l)a;} 

- sin{(a + 2)61 + {(3- l)x} - sin{(a - 2)61 + (/3 + l)a;}. (3. 2) 

This is easily verified, and there is a similar identity with cos(a.9 + /3a;) on the left and all cosines on the 
right. 

Next we require the following integrals over [0,oo) with 9 as in (|l. l| ), F E i(0,oo) and F{oo) — 0: 
I{F,a,(3) = j F{x)sm{a9 + l3x)dx 
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J{F,a,(3) = j F{x)cos{a9 + (3x)dx 
L{F,a,P) — / F{x) cosxain^ 6 coa{a9 + Px)dx. 



It follows from ( [l. 5| ), l|) and an integration by parts that 

{as + l3)I{F,a,f3) 
= jFsHaO + Px)iaO'+P + as-^^coaxsin^O)dx 
= F{0) + J{F\a,P) + as-'^K{F^,a,(3) 

and similarly 

{as + P)J{F, a, 13) = -I{f' , a, ^) + as-'^L{F£,, a, (3). 



Finally in these introductory formulae, it follows from (3. 2) and (3. 3) that 



8K{F,a,f3) = 2/(F,a,/3+l) + 2/(F,a,/3-l) 

- I{F,a + 2,f3+l)- I{F,a-2,(3^l) 

- I{F,a + 2,(3-l) - I(F,a-2,(3+l) 



(3. 3) 



(3. 4) 



(3. 5) 



(3. 6) 



8L{F,a,(3) = 2J{F,a,P+l) + 2J{F,a,P-l) 

- J{F,a + 2,f3+l)- J{F,a-2,f3-l) 

- J{F,a + 2, f3 - 1) - J{F,a- 2,13+1). 



(3. 7) 



We can now return to the integral in ( |1. which we denote by Iq. Then, by (3. 1), we can write 



2/o-/(e,2,l)+/(C,2,-l). 



(3. 8) 
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By ( |3. 4| ) and ( |3. 6| ), the two /—integrals here can be expressed in terms of / and J integrals with 



integrands containing ^ and These last integrals converge more rapidly than those in (3. 8) for cases 
such as 

^{x) = (const.)(l + 2;)-'^. (3. 9) 

Repetition of the algorithm 4| )-( |3. 7| ) accelerates the convergence by introducing integrands with 
higher derivatives and higher powers of ^. 

Certain values of fi have to be excluded to avoid a zero factor as + f3 on the left-hand side of ( |3. 4 ) 
and (|~|). Thus, with a = 2 and /3 = -1 in ( |3. Sj ), we exclude 



(3. 10) 



at the first implementation of the algorithm. At the next application of the formulae ( 3. 4| ) and (3. 5) 
the value 

^ = 1 (3. 11) 



is also excluded and, at the third application of (3. 4) and ( p. 5 ), the additional values 



1/36, 1/16, 9/16, 9/4. 



(3. 12) 



We note that there is an overlap of these values and the iV^/4 values discussed at the end of section 2 
for the validity of ( [Til ) when M > 1. 
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4 Implementing the algorithm 

In this section we show how the integration algorithm from section 3 is used to identify points of spectral 
concentration for potentials of the form 

-c(l + x)"" cosx (4-1) 



given by ( |3. l| ) and ( 3. 9| ), where c > and a > 0. The algorithm is implemented using both symbolic 



methods and numerical approximations. 



We recall from (3. S) that the integral which appears in the formula for p is the sum of two integrals 
/(^, 2, 1) and /(^, 2,-1). The first part of the algorithm consists of a procedure to improve the convergence 
of these integrals. We focus first on the case a > 2 and then comment on the procedure that we have 
been forced to adopt for smaller values of a. 

First an acceptable order of convergence e is decided upon: we have chosen e — x^^ . A purpose 



written Mathematica code is used to repeatedly apply the integration by parts formulae ( 3. 4 )-( 3. 7 ) to 



the integrals whose integrands are larger then e. Starting with the integrals /(^, 2, 1) and /(^, 2, —1), the 



integration by parts formulae generate integrals of the types /, J, K, L (cf. ( 3. 3 )) whose integrands have 
smaller order than the integrand in ij) together with terms that do not involve the variable x. This 
procedure is repeated until all integrands have order less than or equal to e. Thus Mathematica is used 
to generate a symbolic formula which consists of terms that do not depend upon x, denoted by C{fj.), 



together with a sum of integrals of type ( 3. 5 ). Next a Mathematica code is written to parse the formula 
and reconstruct the integrands. The symbolic formula is finally converted into a Fortran 77 function 
which for convenience we denote by F{x, p). 

The next task is to evaluate the integral ( |1. 1| ) numerically from the improved integrand for each value 
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of /i under consideration. As the integral Iq depends upon 9, which itself is a solution of the differential 



equation (1. 5), this is done by solving the system 



^ T( ^\ ( 



(4. 2) 



over [0, X] for some large X, where Io{x,fi) is the integral over (0,a;) in (|l . l| ) , subject to the initial 
conditions 

/o(0,/^) -C(/i), e(0,Ai) = 0, 

typically X=100 has been used in our example. 

There are several practical difficulties encountered in performing the above tasks. First the improve- 
ments in the convergence of the integrals brought about by the integration by parts generate a large 
number of integrals /, J, K, L as indicated by the right-hand sides of ( |3. 6|) and ( |3. 7|) . Numerical 
inaccuracies do not allow us at the moment to improve the integrands beyond a certain point. The 
inaccuracies arise from the K and L type integrals when we attempt to improve the integrand beyond 

The combined integrand F(x,fi) therefore involves 



//// /// 



(4. 3) 



and the excluded values of fj, are ^ and 1 as stated in (|3. 10| ) and ( |3. 111) . We have also tried to perform 



the numerical integration by extending the system of differential equations ( 4. 2 ) to one in which each 
entry is only one of the integrands of /. J, K, L together with the defining equation for 9. However 
this procedure, which has required a more sophisticated parsing routine to be written, has produced no 
significant improvement in the results. 



When a > 2, the order of convergence e — x ^ is achieved by the indicated by (4. 2). When 
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a < 2, we continue to use the same F(x,fi) but with a consequent increase in the truncation error. 
However, the increased error leads to two comphcations which become more serious as a decreases. The 
first comphcation is that ( |1. 1| ) only gives the approximate location of the spectral concentration points. 
The second and more serious complication is that spurious maxima of p are produced, the more so as a 
decreases, and therefore it is necessary to identify the true maxima. We resolve these difficulties by using 
a transitional property of 0{x,fi) which we now describe. 

As reported in section 2] ( see also Q), spectral concentration at a point /io is indicated by a 
certain transitional behaviour of 9{x,p) as /i increases through /ip- Let {xi,X2) be an interval in which 
the Sturm-Liouville coefficient ^ — q{x) < and fi and fi be suitably close to a point of fiQ of spectral 
concentration with /i < < ^ . Then we expect the integral in ( |l. l| ) to be large and negative- and 
therefore producing spectral concentration- if 

{N+]^)^l<e{x,^io)<{N + l)^ (4.4) 

in {xi,X2), where A^(> 0) is an integer. We find that, particularly in situations of sharp concentration, 



( 4. 4 ) is realised with the following features. 

1. 9{x, fi ) and /i ) are close together for < a; < a;i with their values at Xi close to [N + i)7r. 

2. 6{x, p, ) — 9{x, ^ ) is close to tt for x to the right of 2:2, with /i ) close to {N + |)7r. 

3. 9{x, fio) is close to {N + l)7r. 

Thus the graph of 9{x,p) undergoes a rapid transition in {xi,X2) as /i increases from /i to /i . The 
transitional behaviour is illustrated by the graphs in the next section. 

We return now to F{x,fi) in ( [4. 3| ) and the case a < 2. The first complication, concerning the 
approximate location of spectral concentration points, is resolved by using this approximate location as 
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the starting point of a search range for /i within which ( |1. 5[ ) is solved numerically for 6{x, ii). Then the 
points fiQ at which the transition occurs can be located more precisely. We can also use the transition 
property of 6{x,fi) even when a > 2 to verify independently that true spectral concentration points 
obtained precisely as the local maxima of p are not artifacts of our methods. At the same time, any 
apparent maxima which are not associated with the transition property are rejected as spurious. We find 
that, as a decreases, the number of spurious maxima increases significantly making the task of rejection 
a major (but unseen) part of our work. 

Finally in these comments, we note that there is a further complication as a result of the excluded 



values ^ = in ( |3. lOD -( prT^ ). Although is is clear from the integration by parts formulae 

( [3. 4| )-( |3. 7| ) that our methods of evaluating p must exclude these points, the numerical realisation 
of our algorithm produces unreliable results in neighbourhoods of these points. Again we test these 
neighbourhoods for the appearance of the transitional behaviour of 0{x,p). We refer also to ^, section 
6] for a similar use of 6{x,fi) in a different but related situation where a direct formula for p poses 
difficulties. 



In the case of ( [f. the intervals {xi,X2) in ( 4. 4 ) are approximately 



((2r + i)7r, (2r + ^)^) (r = 0, 1, 2, ...). (4. 5) 



We then denote in 

Ai(c, N) [r = 0), iy{c, N) (r - 1), C(c, N) [r - 2) 

and, in the next section, we record our findings concerning the location of these points when c(> 0) is 
regarded as a parameter. We shall also comment on higher values of r as appropriate. 
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5 Examples 



5.1 Example : a = 2 in ( [4. 1|) 



Here q(x) is L(0, oo) and ( |1. 1[ ) is valid for all /i > 0. We apply the algorithm as described in ( 4. 3| ) and 



we have to exclude values of /i near to i and 1 in the computation of p (/i) . In Table 1 we list the spectral 
concentration points which we have located as giving local maxima of p (/i) except that, as mentioned 
in section 4, we have located points near to j and 1 by identifying the value of p which is associated 
with the transitional property of 9{x,p). Thus, for example. Figure 1 gives the graphs of 9{x,p) which 
identify /^(49.26, 1) = 0.25 and provide part of the evolution of p{c, 1) as c varies which is summarised in 
the second column of Table 1. We point out that the graphs represent 9 (modTr), whence the repeated 
cut-offs at the ordinate tt. 

When c has the particular value 122.1, we find that 

^(c,2) = i/(c,3) = 0.51. 

Thus the two local maxima of p arising from /i(c, 2) and i/(c, 3) coalesce or, equivalently, two different 
intervals {xi,X2) ( with r = and r = 1) simultaneously make relatively large contributions to the 
integration in ( p.. 1[ ). One result of this coalescence is that v{c,'i) becomes v{c,4:) when c > 122.1. Thus, 
the last-line entry 0.49 in Table 1 is in fact j^(125,4). 
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c 


Ai(c,0) 


m(c, 1) 


/i(c,2) 


m(c,3) 


iy{c, 3) 


2 


0.45 










7 


0.26 


2.05 








9 


0.08 


2.41 








20 




2.10 








30 




1.87 


5.28 






40 




1.29 


5.51 






50 




0.15 


5.10 




0.68 


60 






4.81 




0.72 


70 






4.73 


9.75 


0.70 


80 






4.45 


9.!).") 


0.67 


100 






3.11 


9.26 


0.60 


120 






0.80 


8.9 


0.52 


125 






0.10 


8.84 


0.49 



Table 1: a = 2 
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5.2 Example: a = 1 in (4.1) 

Now M = 1 at the end of section 2, and ( |l . l| ) is valid for all /i > except /i ~ 1/4. We have followed 
the procedure for a < 2 as described in section 4, in which the excluded values are /i = 1/4 (again) and 
/i = 1. We find that spectral concentration points occur in greater profusion than for a = 2, and our 
results for c < 100 are given in Tables 2 and 3. Coalescing points occur as follows: 

/i(c, 1) = iy{c, 3) = 1.98, (c = 24.25) 
Ai(c, 1) = ,^(c,4) = 0.53, (c = 31.2) 
H{c, 2) = zy(c, 5) = 4.02, (c = 64.6) 



as well as 

/i(c, 3) = i/(c, 8) 10.06, (c = 106). 

As c increases through these respective values, the oscillation number N for i'{c,N) or £^{c,N) increases 
by one, because of the transition tt in the values of caused by the fj, spectral concentration points. This 
increase is indicated at the tops of columns in Table 3. Also, in Figure 2, we give the 0— graphs which 
show the transitional behaviour twice to illustrate the coalescence of fi{c, 2) and ^{c, 5) when c = 64.6. 
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c 


A^(c,0) 


A*(c, 1) 


m(c,2) 


m(c,3) 


/i(c,4) 


2 


0.5 










3 


0.48 


3.05 








4 


0.41 


2.73 








5 


0.26 


2.68 








6 


0.03 


2.68 








10 




2.76 


6.77 






15 




2.97 


6.81 






20 




2.61 


6.84 






25 




1.85 


7.09 


12.51 




30 




0.80 


7.36 


12.73 




40 






7.35 


12.82 


18.38 


50 






().43 


l:-!.3(i 


20.41 


60 






4.89 


13.90 


20.55 


70 






2.89 


13.92 


20.79 


80 






1.75 


13.36 


21.39 


90 






1.55 


12.35 


22.06 


100 






1.31 


11.01 


22.47 



Table 2: a = 1 
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c 


u{c,2) 


zy(c,3/4) 


z.(c,5/6) 




^(c,4/5) 


2 












3 












4 












5 


0.70 










6 


0.75 










10 


0.72 










15 


0.59 










20 


0.41 


2.22 






0.69 


25 


0.14 


2.06 






0.64 


30 




2.05 


4.05 




0.57 


40 




1.83 


4.10 




0.34 


■50 




1.40 


4.15 




0.19 


60 




1.00 


4.08 




0.04 


70 




0.43 


3.88 






80 






3.60 






90 






3.24 






100 






2.82 


10.14 





Table 3: a = 1 
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'mu=3.97' 

'mu=4.02' 




2 4 6 8 10 12 14 

Figure 2: 6 graphs for /i(64.6, 2) = iy(64.6, 5) = 4.02 
5.3 Example : a = | in (4.1) 

We have intimated in sections 2 and 4 that slow decay spawns theoretical and computational complexities, 
the latter including the task of segregating the large number of spurious and actual local maxima of p . 
This task is expensive in computer time and therefore we have restricted the range of c in this part of 
our investigation to < c < 30. Our findings for ji, v and ^ points are summarised in Table 4, in which 
figures are given to more than two decimal places when necessary to distinguish between closely situated 
points. 

In particular, we find that there are two values of c and \x in whose neighbourhoods several very close 
(but apparently not coalescing) points of spectral concentration exist. The values are 

c=3, /x = 0.49; c = 5.5, /x = 0.66. 
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In the case of c = 3, for example, we have evidence for the existence of at least seven such spectral 
concentration points, corresponding to < < 6 in ( [4. 5| ). In Figure 3, we exhibit this evidence for 
< r < 4. The figure shows 0— graphs for a decreasing sequence of values of /i (all close to 0.49) where 
the transition occurs for r = 0,1,2,3,4 in turn. We make a further comment on this matter in section 
6.3 below. 

For c in the stated range (0, 30) we have identified one pair of coalescing points which, as in section 5.2, 
are associated with a change in the oscillation number N, as follows: ^{c, 1) — i/(c, 3) = 0.67 (c — 24.5). 
Further, in addition to what is reported in Table 4, we have found that, for 15 < c < 24.5, there is a 
v{c,A) which is almost identical to /i(c, 1). Thus, for c in this range, we have a spectral concentration 



point which enjoys contributions from both 7' = and r = 1 in (4. 5). This linkage between and v 
points is a new feature for a = 1/2 and it poses the question whether there is a theoretical explanation. 



c 






Mc,2) 


u{c, 1) 


u{c,2) 


i/(c,3\4) 


i{c,2) 


?(c,3) 


€(c,4) 


?{c,5) 


2 


0.61 






0.56 














3 


0.50 






0.487 


0.73 




0.4869 


0.62 






4 


0.26 








0.72 






0.66 


0.78 




6 










0.575 


1.83 




0.569 


0.71 




8 




3.49 






0.33 


2.08 




0.331 


0.58 


1.85 


10 




3.61 






0.09 


2.11 








2.01 


15 




3.11 








1.842 








1.841 


20 




2.0 


8.82 






1.2986 








1.2984 


25 




0.51 


9.02 






0.59 










30 






8.72 

















Table 4: a =1/2 
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1 


1 1 


1 1 






1 

'mu=0.51' 


— 












'mu=0.4885' 














'mu=0.487' 














'mu=0.48685' 














'mu=0.48681' 
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Figure 3: 9 graphs for a = i, c = 3 



6 Concluding remarks 

6.1 Absence of spectral concentration 

In section 5], we made the following conjecture concerning the appearance of spectral concentration 
for (1. 2) and ( |1. 3| ) when q{x) has the form cQ{x) (c > 0) such as we have in ( p~^ ), ( p~8| ) and (g^). 



Conjecture 6.1 Let Q{x) < in some interval (0, xq)- ^ei change sign one or more times as x 

increases with Q{x) finally decaying to zero as x oo and Q{x) e L(0, cx)). Then there is a number 
Co (> 0) such that spectral concentration does not appear for any /i > when c lies in the range (0, Cq). 



For the example ( |l. 7j ) we showed in Q that the conjecture is true with 0.28 < cq < 0.29 and, more 
generally, an affirmative answer to the conjecture has been given recently for the special case of ( |2. l| ) in 
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which x^{x) is L(0,oo). Thus is covered provided that a > 2. The conjecture remains undecided 

when 1 < a < 2 in ( [4. although our computational findings for a = 2 indicate that cq exists with 
1.0 < Co < 1.1. 

6.2 More general potentials 

There is no difficulty in principle in extending our methods for ( |2. 1| ) to potentials which are a finite sum 

L 

q{x) = ^^dx)piix) (6.1) 
1 

and the pi{x) have different periods lui, provided that the uji are mutually commensurable. The pi{x) 



would then all have a common basic period lu, and the set cr in (2. 6) is modified to 



cr = {iVV/tj2;iV = l,2,...}. (6. 2) 

However, in the absence of a special form such as ( |2. l[) or ( |6. l| ), it is no longer clear what can be said 
in general about the absolutely continuous nature of p when q fails to be L(0, oo). There is no simple 
exceptional set a such as ( |6. 2)) because in [|l^, Section 4.4], ||ll[ and |2^ it is shown that, given any set of 
isolated positive real numbers /i„ and given a (0 < a < 1), there is a potential q{x) — 0{x~'^) (x — > oo), 
such that the spectral function has discontinuities at the /x„. 

6.3 Higher values of r in (4.5) 



We have focused on r = 0, 1, 2 in (4. 5) in order to establish our computational methods, and consequently 
we have located mainly /u, v and ^ points of spectral concentration. Thus we have largely confined our 
search for the transitional behaviour of 9{x,ii) to the a:;— range (0, 18). There remains the question, on 
which we have touched in the remarks relating to Figure 3, whether 6{x, p.) possesses the transitional 
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property also for higher values of r and, more particularly, whether an infinity of values of r is involved 
in this way for some fixed value of c. Thus we have the theoretical question of whether an infinite set 
(bounded or unbounded) of spectral concentration points can exist for some c. All that is known is that 
an unbounded set cannot occur when a > 2 because it is shown in H that the set of spectral concentration 
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